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Abstract 

We present analytical methods for investigating the interaction of two heavy 

quarks in QCD3 using the effective action approach. Our findings result 

in explicit expressions for the static potentials in QCD3 for long and short 

distances. With regard to confinement, our conclusion reflects many features 

found in the more realistic world of QCD4. 
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I. INTRODUCTION 



Gauge theory models in other than 3 + 1 space-time dimensions have been a center of 
interest for many years. In the present article we want to explore QCD 3 and work out its 
similarities and differences in comparison to QCD4, the presumably correct theory of strong 
interaction. Our investigation follows in many respects the work of Adler and Piran ||T| in 
QCD4. Their starting point, the renormalization group approach H, is, however, bound to 
fail; there is simply no renormalization group in QCD3 nor is there a Callan-Symanzik (3- 
function, etc., well-known attributes of QCD 4 ; our theory, QCD 3 , is super-renormalizable. 
However, the infrared problem, which we are about to analyze, is shared by both theories. 
Now, since the incompletely understood gluonic vacuum structure is at the heart of both 
QCD 4 and QCD 3 , we must rely on some more or less reasonable effective action models. In 
Adler's case it is the leading-log model while ours might be termed the leading-root model. 
The latter will be represented by the one-loop effective action with constant color magnetic 
background field. Among many interesting features revealed in our work is the QCD3 
vacuum acting like a dielectric medium, the elliptic shape of the confinement region and 
its scaling properties and, finally, the behavior of the static potential between two massive 
color test charges for large and small distances. In Section II we present the essentials for 
the calculation of the one-loop effective action. In Section III we focus on the large distance 
(confinement) problem and give an expression for the linearly rising potential plus correction 
terms. In Section IV we treat the short distance domain and derive the classical formula 
for the interaction of two static charges augmented, again, by correction terms. Section V 
summarizes our findings. 



II. COMPUTATION OF THE ONE LOOP EFFECTIVE LAGRANGIAN 

Since much work has been invested in the calculation of effective Lagrangians, we will 
use some short cuts to quickly reach our present goal. The Lagrangian for a pure SU(N) 
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gauge field theory is 



(1) 



where F afiu = d^A au - d v A a " + g f abc A h ^A c 



(2) 



Working in d = 3 dimensions we list the following dimensions of various quantities, some of 
them to be introduced at a later stage: [g] = [Q] —[A ] = m 1//2 , [J°] = m 5//2 . 
Next, the gauge field is decomposed into 



(3) 



where b atl represents the fluctuating Yang-Mills field. The external field configuration A^f 
that probes our QCD vacuum is the commonly chosen color magnetic background field 



^0 ^ 



_ _1 pM^™ rp 



B 
-B 



(4) 



where F^ v = Fg u V a , A% = A»V a , (V a ) 2 = 1. 



(5) 



The magnetic field points along a fixed unit direction V a in color space. Inserting the 
parametrization (|3|) into the Lagrangian ([!]) we obtain 



- (DtD%) + g f abc F£\ \T + 0(b :i ), 



abc rpb 



(6) 



where D$ = (d^6 ab + g f acb A 



B/J 



The Lagrangian (KJ) together with the background gauge fixing term, 



C GF = ^(D^DlW 



(7) 



(8) 



connects us with the effective Lagrangian via the functional integral 



iV'exp 



i / d 3 x£ eff 

N [ £>6A FP exp 



i / d 3 a; (£ + £ 



GF, 
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Using the Feynman gauge, a = 1, we obtain 



£ + £ 



GF 
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To compute the two terms in the square brackets of (|10|) we follow the author of ref. |§ and 
so reproduce the result (pay attention to signs!) 



£ + £ GF = -\B 2 + J2 Ib^g^-d^b^ 

T b £C v 

+ £ W£*(D 2 g^ + 2igQ a F, v )W^ 



(11) 



where = d» + \gQ a A^ = - \gQ a F» v x v 



andCy = {T b \ [V,T b ] = 0} , V = V a T a , 



(12) 



where the T a 's denote the standart SU(N)-generators: [T a ,T b ] = if abc T c , Ti{T a T b ) 
\5 ab . The remaining generators (jL Cy can be expressed in terms of eigenvectors of V with 
eigenvalues Q a . Their fluctuating Yang-Mills field components b atl form certain complex 
linear combinations W£ depending on the choice of V in color space. 

At this point we perform a Wick-rotation so that all subsequent calculations are done 
in Euclidean space: 



iV'exp 



d 3 x £ E ff 



N" J Vb A FP exp J d 3 x (£ E + £ EGF 
N'" J VW*VW A FP exp J d 3 x(-\L 

+ E W*(D\„ + 2igQ a F flv )W ( 



(13) 



Due to the fact that the fields W* and W have to be treated independently the Faddeev- 
Popov Determinant appears in the form 

A FP = A^A^p = Bet(-D 2 )Det(-D 2 ) 

= exp[21nDet(-D 2 )] = exp[2Tr ln(-D 2 )] 



exp 



J d 3 x2tr ln(-D 2 ) 



Now we can write (1131) as 



N' exp 



d 3 xCf 



N'" exp 



d 3 x(-±5 2 + 2trln(-D 2 )) 
x f VW*VW exp f d 3 x J]W^pV 

J L a 

+2i^Q a F Mi/ )W^] . 



Here we meet the standard functional integral 



/ 2W*2W exp - J <^xW;M^W v 
= exp[-Tr lnM Mi/ ]. 



(14) 



Hence eq. fll4]) takes the form 



N' exp 



d 3 x Cf 



N'" exp 



/" d 3 x(-±B 2 + 2trln(-Z> 2 



-trln(-D 2 ^-2i^F^)) 



from which we read off the expression for £| ff : 



eff _ 1 D 2 



5 2 + 2trln(- J D 2 ) - ti\n(-D 2 S^ u - 2igQ a F f 



/IV J 



(15) 



(16) 



There are several ways to compute the various traces in this expression. We prefer the 
^-function regularization method 0. The result is 
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hB 2 + \gB\' 



(V2- 



2?r 



47T 



■« 



1 3/2 



(17) 



In ref. || and we found a proper-time calculation of expression (|T6|). We agree with the 
result contained in ref. |3[]. The numerical value of is ~2. 61238 so that the sign in front 
of \gB\ 3 / 2 is indeed positive. Hence C eS takes a maximum (V eS a minimum) at a nonzero 
value of the background field: 

(V2-l^ 



\B l 1/2 



3/2 A 

y 4,,- 



47T 



£IQ«I 2 - 



(18) 



The gauge invariant generalization of ([T7[) can be obtained by the replacement £> 2 

^«)2 _ (£a)2 = ^. 



1 r^ajiv rpa 

2 \w 



( 



c 



cff 



1 



(19) 



where k 1/2 = |5 CX | 1/2 = eq.(18). 



(20) 



So we obtain a gauge field condensate T = k 2 due to radiative corrections, just like in four 
dimensions, which determines the interesting features of the model. If we choose SU(N=3) 
as our gauge group and the unit color vector pointing along the three-direction we find 
k 1 ' 2 = 0.37245... g 3 / 2 . 

At this stage it is important to point out that we are only dealing with the lowest (first) 
order loop approximation. Higher order loop calculations will certainly modify the position 
and the shape of the extremum, i.e., the value of k. However, concordant with our own 
approach there are other strong indications for gauge field condensation as implied, e.g., by 
the so called average action method advocated by the authors of ref. Thus, we expect 
our leading-root model to describe the exact QCD 3 -vacuum structure at least qualitatively 
accurate. 

Furthermore, we assume that the form of ( |T9"D also holds for static fields which are slowly 
varying in space. 



III. FLUX CONFINEMENT AND THE HEAVY QUARK STATIC POTENTIAL 



In this section we study the statics of two massive test charges at large distances. The 
approximation in which leading QCD 3 radiative corrections are retained is given by the 
effective Lagrangian (|H|). Following Adler we write for the potential of static, infinitely 
massive test charges [J® = ( Jq , 0)] 



Ktatic = -extr A a | J d 2 x (C eS (AD - A a J£)} 



(21) 



Limiting ourselves to the case of quark-antiquark source charges we set (R = 2a) 



J* = Qq a [5 2 (r - ax) - 5 2 (r + ax)} 



(22) 



where q a is the internal unit vector in color space. Similarly we introduce scalar and vector 
potentials ip and A by writing 



A a = q a V , A"- — q a Aj. 



(23) 



Our variational problem can now be restated in the form 



-extr A , v | J d 2 x (C^F) - ipJ ) 



where 



(24) 



T ' = E 2 -B 2 , E = -V<p , B = e^Aj 



(25) 



With £ eff (jF) given by ([19]) we write more explicitly 



Static = -extrA,^ 

/ / 



d 2 x 
: extr A ,<p 



\ 



1 - 



V 



T_ 



d 2 X -^stat 



(26) 
(27) 



Given £ sta t we now can apply the Euler-Lagrange equations which imply the field equations 
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V-D = J , eijdiEj = = eijdiH , (28) 



T_ 



(29) 



where D = eE , H = eB , e = 1 
The source-free equation (^) for the magnetic field can be satisfied by 

eB 2 = 0. (30) 
Here we have to distinguish three cases: 

(Ia)fl = , E 2 > k 2 

(Ib)fl = , -E 2 < /t 2 (31) 
(II) e =0 , B 2 = E 2 -k 2 . 

For short distances we expect Coulomb-like field configurations with E large and B vanish- 
ing. Hence there should exist a finite region containing the source charges for which (la) is 
satisfied. In this domain we have reduced our original variational problem to a problem in 
nonlinear electrostatics with a field strength dependent dielectric constant: 

V • D = J , CifrEj = 0, (32) 
D — e{E) E , e(E) = l-^, E — \E\ e R. (33) 

Now, in Adler's leading-log model of QCD 4 , it proved very effective to work with a manifestly 
flux conserving quantity. So let us likewise parametrize D by introducing a scalar flux 
function D = /($)• Without going into all the details (and subtleties) of how to arrive at 
the explicit relation between D and $, we just state the result: 

(The authors of ref. missed the important factor |.) The boundary conditions imposed 
on the flux function <&(x,y) are 
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$(|x| < a,y -> 0) = Q, 

$(x > a,y -> 0) = 0, (35) 
$ — > for x 2 + y 2 — > oo. 

For future calculations it is useful to derive some relations between the fields E and D. From 
eqs. ( |3~3"D and still considering the branch (la), where E > n, we obtain 



D=|.D| = |e||.E| = e£;, e= A /l--|>0. (36) 



This equation implies -D = —^[kE x I 2 + or, setting e := 



e 1 t = ^±^+D. (37) 



2 tZ , i Kp 1 



We need ef 9 = - ± W h kD 

1,2 2 V 4 

and select the positive sign to guarantee a single-valued potential. So we have the relation 



E = E(D) = ^ + ^ + kD + D. (38) 

At last we turn to the solution of eq. (P^)- In order to obtain more insight into the behavior 
of the solution, we begin by rewriting the field equation for <&(x, y) in its characteristic form. 
So let us start with 

eijdiEj = 

or = d x E y - d y E x = d x - d y (—} 

Here we employ eq. (|37D and introduce the flux function via relation (|34|) . The result of a 
rather lengthy chain of partial derivatives is given by the exact field equation 

dl + dl + (a- 1)«9 2 ) $ = 0, (40) 



where d n = h • V , ft = 



2 /_2 



v i r +K - D 

and a = . (41) 

2D D 

Letting d\ be the tangential derivative, we can replace the coordinate derivatives (d x ,d y ) in 
©by 

+ dy = df + + first derivative terms 

and so the field equation ([40D takes the form 

(# + a^)$ + 0(3A) = O. (42) 

Now, for weak fields, D — > 0, i.e., away from the charges, we have 

f-r2 — + J D 

/ 2 



lim a = lim , = 0, (43) 

d->o D-*oz + D + y /e» +KD 

and for strong fields, D — ► oo, i.e., close to the charges, we obtain 



lim a = 1. (44) 

The flux equation (flUD is quite similar to the one found in QCD 4 ; it is of degenerating elliptic 
type and has a real characteristic at a surface of constant $, where V$ = 0. Using the same 
arguments as in Adler's leading-log model, we have here the first indication of confinement 
in QCD 3 . Next we want to show quantitatively that, in fact, the total flux between two 
massive color charges is confined to a domain with a characteristic as boundary on which 
$ vanishes. To do so it is useful to reformulate our problem in still another form for the 
equation for For this reason let us go back to eq. fl32|) . We also know that E depends on 
D in a way stated in (|38|): 



K / k 2 



f(D) := E(D) = - + J- + kD + D. 
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Then we obtain 



= d x E y - d y E x = d x - d y (j^f(D) 

= d x (^f(D))-d y (^f{D)). 



Recalling relation ( |34] ) we get 

d x 

+ dy 



d x $f{D) 



_((d x <s>y + (d y m i/2 \ 
((^$) 2 + (d y $yy /2 



(45) 



where ((d x $) 2 + (d y <5>) 2 ) 1/2 = 2D 



After performing the various partial derivatives we end up with the following, still exact, 
differential equation for 







D 



+ 



9(D) 



1 (d x $) 2 + (dy*)'' 



9> 



D d X $dy<S>d Xy <S>, 



where g(D) := 421 - D. 



rip) 



(46) 



(47) 



At this stage we make contact with calculations contained in ref . || . Needless to say, a solu- 
tion of fl4"B] ) is not easily available. However, being interested in the far-field approximation, 
we now study the limiting case of weak fields: 



f(D) = % + D + fa + KD 

k k ( 2D D 2 ' 

-2 +D+ 2( 1 + -- 2 ^ + 0i3) 



so that to first order in D: f(D) = k + 2D. From here we obtain for fl4"T|) 



9(D) 



(48) 



meaning g(D) ~ £>(L>°) 
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Using the value ([48]) in eq. ([46]) we arrive at 



+ 



<9 2 $ 

v 



-d x <$>d y $d xy <$>. 



(49) 



Following the strategies in QCD 4 it is convenient to rescale x and y in terms of dimensionless 
parameters: 



x = Rx , y = R a y. 



(50) 



The authors of ref. M supply arguments as to why the transverse coordinate scales with 



a 



Now, we try the following ansatz: 



$ = $(0) + + J_$(2) + _ _ _ 

R R 2 



Earlier we found {2D) = ((^.$) 2 + {d y <5>) 2 ) l/2 , so that 



1 
2D 



R 2 



+ 



1 



^/3-(^( $(0) + i* 

R2/3 (^$(°)) 2 



(1) 



+ 



-l/3\ 



with <9g$ (0) < 0. 



(51) 



(52) 
(53) 



Again, we skip the details of the calculation for the various partial derivatives in ([491). Then 
we obtain the following differential equation for the flux function in zeroth order (x, y = x,y): 



0= K^Wj^Wj+Kjll^^)^ 

- (d 2 ^)(dy^) 2 - 2K(d X ^)(d X y^). 



(54) 



Another useful form of this equation can be obtained by multiplication with (d y &^) 2 : 



o = a, 



2 \d y &°\ 



(55) 
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We do not expect eqs. (j54|) or (|55|) to be soluble for y) by the separation of variables 

method since there is still the boundary condition (|35"D to be taken into account. Yet, as in 
QCD4, there is hope for the existence of a separable solution y(x, $). Hence our next step is 
to rewrite ([54]) into a differential equation for y. Here are the rules governing how to achieve 
this ($ = 

% = V*\ = %v = - V -f (56) 

eft yg*y| , n Vx*yx 

^XX o ~T A „ 

y% y% 

This set of partial derivatives enables us to rewrite (|5lD in the form 



^r-K-jr = °- ( 57 ) 



Here we can try the ansatz 



y(x,$)=X(x)F($) (58) 



and so obtain for (1571) instead 



c, (59) 



FF' 3 

kX"X 2 = c, (60) 



where c is the separation constant. Integrating (|59| ) yields 

* = -^ + A 'xM + ^ (61) 

With the aid of ( [351) we find for fc 2 : 

<% = 0) = Q : fc 2 = g. 

Since it is easier to work with X(x) in the numerator of (|6lD , we redefine X — > -4 and so 
obtain 
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<t> = Q{l-a iy X{x) + a 3 y 3 X 3 {x)) (62) 

k\ 1 c 

with new constants ai = — — , a 3 = — — . (63) 

Q 6Q 

Now recall that (i) we have to satisfy d y Q < in the confining domain and (ii) we also want 
$ to approach the boundary $ = continuously: d y <§>\ yh = 0. The first condition tells us 
that ai > while the second implies 03 > 0. These considerations lead to the following two 
equations: 

= ®(x,y h (x)) 

= Q(l - ai y h (x)X(x) + a 3 y 3 (x)X 3 (x)), (64) 
= d $1 

= Q(-aiX(x) + 3a 3 ylX 3 (x)). (65) 
The last equation yields the explicit expression 



yhix) = ± ^itxTx)- (66) 



The two signs reflect the symmetry with respect to the x-axis. Substituting (pq) into (| 
we obtain 

4 



~ til 

27 1 

Thus, there is only one free parameter left. So far we have 



(67) 



y) = Q{1 - \byX{x) + ±b 3 y 3 X 3 (x)) (68) 

where b = |ai , or a 3 = |6 3 . (70) 

What remains is an explicit solution for X(x). Because of our redefinition X — > the 
equation following from (|60|) is 

«(l)"x- 2 = c. (71) 
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Introducing X' = ^fp we can cast ([71]) into the form 



dX X k 



(72) 



which can be solved by 



p(x) 



6Qb 3 



X 4 



X-Xn + 



np 



(73) 



' 6QX$b\ 

where X = X(x = 0) and po = p(X Q ). Among the three free parameters b,X ,p we find, 



via equation (|69|), 



1 



yb{x)\x=o 



bX n 



(74) 



that b and Xq merely rescale the size of the confinement domain. Hence solution (|T3|) is 
essentially determined by p . Because of 



Po = p(X ) = p{X{x = 0)) 



'dxY 



2 = 



we see that p is (a) related to the tangent of $ in x-direction at x — 0, 

d x $\ x=0 = Q(-\byX' + pyX 2 X% =0 

and (b) related to the shape of the confinement boundary at x = 0: 

1 



(75) 



d x y h (x)\x=o 



bX- 



:X' 



x=0 



(76) 



Now there are three different cases to be distinguished: 

6QX 5 b 3 
1 Poi = 

K 

(2) p 2 = 

(3) P03 Poi , P02 , arbitrary . 



(77) 



Case (1) has been treated in the literature and has the advantage of being analytically 
soluble. We regard this solution as unphysical, since, looking at (|T6j), yb(x) does not behave 
smoothly at x = 0; this solution never yields an extremum with regard to, e.g., the energy 
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density, assuming an underlying variational principle. On the other hand, case (2) contains 
a physical, smooth boundary, d x y^\ x= o = 0, and is our preferred solution. Its disadvantage 
is that it cannot be solved analytically. Case (3) is neither analytically soluble nor is it 
physical and so will not be considered any further. 

Without going into further details we now present our solution of eqs. (Q,^) for the 
case (1): 

V3 



[a — x 



)2/3 



k y 



27Q (a - x) 2 
3 fCh ]: '' 



2/3 



(78) 
(79) 



(where we have switched back to our original coordinates x, y — > ^x, ^j^y). We have 
assumed x, y > 0; otherwise we would have to use moduli. Note the scaling behavior of y 
with a |-power. 

Finally we come to the calculation of the static potential. This is achieved with the aid 
of the formula 

D 



Static = Jd 2 x J E(D') dD' 




D 

d 2 x J dD' 





d 2 x\ 



K IK 

2 " 



— + kD' + D' 
4 



2_ hf_ 

3^VT 



kD 



K 



12 2' 



(80) 



The surface integration has to be performed over the confinement region defined by ([79 

Ktatic. = / d 2 x 



a-e Vb(x) 

i] = J dx J dy [%) 

-a+e -y b (x) 
' [' < I 





4 J dx J dy 
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Restricting ourselves to weak fields, i.e., to large distances (confinement or infrared domain), 
we approximate the integrand in ( ftfty by 

[%] = kD + D 2 + . . . , 

where D = ^[(d x &) 2 + (c^) 2 ] 1 / 2 in which we substitute the perturbative ansatz (|5T|) and 
rescale according to (|50|). An intermediate step on our way to Viatic is 



static 



h + L 



2 • 



(81) 



where 



a-e Vb{x) 

Ii = —2k J dx J dy 





-2k J d% J dyd y $ 





a— e 

— 2k J dx 





=0 =Q . 



a, — t 

2kQ /" dx 



lim 2«Q(a — e) = kQ-R 



So indeed we have linear confinement! 
The correction term I 2 turns out to be 

12 



6 -JQ 7 o-i/ 



(82) 



(83) 



Here, then, is our final expression [case(l)] for the static potential of two oppositely charged 
massive test sources at large distances: 



Static = KQR + ^^QB R l/3 _ 



/3 



+ 



(84) 



5 v 35 V k 

Note the linear rising with distance as it is familiar from QCD 4 . 

Now we turn to case (2) with po = P02 = 0, which guarantees a smooth boundary, 
X'(x)\ x=0 = 0. The equation to be integrated follows from (|73|): 

dX 



dx 



^ X 2 \/A — X . 



K 



(85) 
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The solution is easily found: 
x = ± 



QQ(bX ) 



X(x) 
Xo 



X(x) 

x 



1 X(x) 1 
h arctan 1 / —rz 1 



X 



(86) 



This equation is transcendental with respect to X(x), i.e., neither X(x) nor Q(x,y) can be 
written in an explicit analytical form. It is, however, possible to find an approximate solution 
which is sufficiently close to an exact solution and which maintains the smooth boundary 



yb(x). Defining t = X{x)/Xq it turns out that y/t — 1/t + arctan \/i — 1 can be excellently 
approximated by arctan y/t q — 1 with q ~ 3.301..., whereby q = 3 already represents a 
fairly good approximation. So we can write instead of (f 



x — ±, 



6Q(&X ) 3 



arctan , 



x 



- 1 



or X(x) = Xq I 1 + tan 



'6Q(6X ) 



■ x 



/v 



1/9 



(87) 
(88) 



With the aid of eq. (H), i.e., employing the condition 



= y h (x = I] 



1 + tan 2 



bX 

we obtain bX 



l 6Q(bX y ! 
k 2 

/ 2\ 1/3 



-1/9 



6g 



Hence our approximate solution for the flux function is given by (x,y — > x/R,y/R 2 ^ 3 ): 



$®(x,y) = Q 



/n 2\ I/ 3 



V 16Q I R 2 ' 3 



(l+tan 2 ff) 



1/9 



/t7T 2 ?/ 3 



, 2 7r a; 



+ i2g^( 1 + tan - 



3/9 



(89) 



with the confining boundary 



2/b(^) = 2 ( — 



,2/3 



(1 + 



tan 



2 7T X 

2 a 



1/9 



(90) 



Again, the scaling behavior x oc a and y oc a 2 / 3 is visible. 



With the aid of the flux function we now turn to the calculation of the static potential. 
According to eqs. ( |3T| , |32"D there is no change in first (linear) order oc kQR; for the calculation 
of this term the explicit form of $ is not needed. The correction term yields 



a-e Vb(x) 

I 2 = - J dx J dy 





K 



4 (<9 y $(°)) 3 



and after performing the various integrals we end up with 

Static = kQR + 2.17.. {f^R 1/3 - 0.36.. {/^ iT ^ 



(91) 



(92) 



IV. SHORT DISTANCE BEHAVIOR 



Let us recall that in four dimensional space the static potential can be evaluated by 



^static = -{extr^W - AVcoulomb) , 



(93) 



where W 



d 3 x 



|6o(V^ln(^)-^ 



(94) 



By means of a cleverly chosen rescaling of coordinates, fields and charges, Adler succeeds in 
splitting up the effective Lagrangian into a classical part and extra terms due to quantum 
corrections containing a dimensionless running coupling ({R). This function goes to zero as 
R approaches zero. Hence limiting oneself to short distance behavior, ip and Viatic can be 
expanded in a perturbation series around the classical Coulomb solution with Q{R) as the 
small parameter. 

In QCD 3 we found for the effective action 



W = J d 2 [C eS (E) - <pj 



(95) 
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Here it is obvious that the classical part oc (Vy?) 2 is already separated from the part contain- 
ing the QCD corrections. Hence, rescaling the various parameters does in no way improve 
the situation. In fact the static potential approaches its classical limit in a natural manner, 
since close to the test charges we have E ^> k. Thus, for short distances, R — > 0, our zeroth 
order approximation is sufficiently descibed by electrostatics: 



Static = -extr v J d 2 x \{Vip) 2 - ipJ ( 



From (|96|) we obtain the Poisson equation, 



(96) 



V 2 y? = —J = —Q[5(r — ax) — 5(r + ax)] , 



(97) 



which has the well-known solution 



Q . \r — ax\ 

— In r 

2tt 



r + ax\ 



Q y/(x - a) 2 + y 2 
27r ^(x + a) 2 + y 2 



The Coulomb counter term becomes 



AV C = extr^jl d 2 x [|(V^) 2 - yJ x 

+ J d 2 x [I(V^) 2 -^J 2 " 
Ji = Q5(r — ax) , J 2 = —Q5(r + ax) 



(98) 



and the solutions of the corresponding field equations are 



y?i = in A \r — ax , ip2 = — InAr + aa; 

2n 2% 



(99) 



where the arbitrary parameter A has dimension [L] 1 . Using the Poisson equation ( j97|) in 
the action we find for W and similarly for AV C 



ir = / d 2 x 

d 2 x 



i(V^) 2 -^J 



■\ I d 2 x(pJ . 



Working with this expression for the action we obtain for Viatic 
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Static = \ J d 2 XLfJ - \ J <\ 2 XifiJi - \ J d 2 XLp 2 J 2 



2tt 



(100) 



This, then, is the leading short distance behavior of the static potential in QCD3. In the 
sequel we will demonstrate that expression (|100|) also shows up when we now compute V^ tat i C 
by means of the formalism we developed for large distances. To do so we have to return 
to the exact quasilinear, second order differential equation ([46]). Expanding the coefficient 
function g(D) in terms of D near the sources where D 3> k, we obtain 



g{D) = \^D + - + 0{D-^) 



(101) 



The relevant ratio g(D)/D is thus of order 0(D 1 ^ 2 ); hence we are permitted to omit it 
for short distances, obtaining, not surprisingly, Laplace's equation as approximation of eq. 



(d 2 x + d 2 )$ = 0. 



(102) 



The well-known solutions satisfying the boundary conditions (IJH) are 



<3>(x, y) — $ 



cl 



Q 



71 



y y 

arctan arctan ■ 



x — a 



x + a 



This expression for our flux function can be used to find 



D — - 

2 



V 



'y 



Q_(ti_ti 

2tt \n r 2/ 



(103) 



(104) 



Again, we come to the hardly surprising conclusion that the classical linearization is sufficient 
for treating the short distance behavior of two static color test sources. 

At last we turn to the calculation of the static potential, which can be achieved with the 
aid of formula (EC 



Static = / d 2 x[\D 2 + f + f D + 0{D^ 2 )} 



_. yD\ , yV'^ , yu , 
■ static 1 static 1 static 1 



(105) 
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Of course, V s ^ tic yields the classical potential and, with due consideration of the Coulomb 
counter terms, reads, as before, 

KuL = ^ In Ait:. (106) 

The other two remaining terms in (|105|) allow us to augment the classical potential by 
correction terms. Needless to say, the respective integrations have to be performed with 
great care ||. The leading order contribution to the classical potential comes from Vj^?, 
while V^ tic provides us with a correction term oc R. Our findings for the potential of two 
static color charges at short distances can be summarized in 

Q 2 

Ktatic = — In Ai? 

271 



2C 

+ — kQR, (107) 

7T 

where ty'(x) denotes the derivative of the psi function with ^'(1/4) ~ 17.1973 and G is 
Catalan's constant, G = 0.915965. . . . With these numbers we can rewrite (|107|) in the final 
form 

Static = 77- In XR + 0.054.. JkQ 3 R 1/2 + 0.583.. kQ R . (108) 

Hence in addition to the dominant classical potential, there exist subdominant contributions 
behaving like i? 1//2 ,i?, . . . vanishing as R approaches zero. Equation ( |108j ) should be read 
side by side with Adler's formula (40) of ref. [ffO . 



V. CONCLUSION 

The results in this paper point to great similarities between QCD3 and QCD4. Despite 
major differences in physical and analytical details we find that in both theories radiative 
corrections to 1-loop order spontaneously generate a gauge field vacuum condensate leading 
directly to a confining theory. We find it interesting that quantum contributions arise 
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from the employment of purely classical differential equations. In this way large distance 
as well as short distance correction terms to the classical potentials were found. Without 
overestimating the importance of low-dimensional field theories, our calculation may lend 
some further insight into the mechanisms of classical approximations of QCD 4 . 
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